Relativistic effects in proton-induced deuteron break-up at 
intermediate energies with forward emission of a fast proton pair 
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Recent data on the reaction pd {pp)n with a fast forward pp pair with very 

o : 

^ ■ small excitation energy is analyzed within a covariant approach based on the Bethe- 



Salpeter formalism. It is demonstrated that the minimum non-relativistic amplitude 
is completely masked by relativistic effects, such as Lorentz boost and the negative- 
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^ ' energy P components in the Sq Bethe-Salpeter amplitude of the pp pair. 
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^ , The investigation of hadronic processes at high energies, such as reactions of protons 

O 

scattering off deuterons, provides a refinement of information about strong interaction at 
short distances. Nowadays, large research programs of experimental studies of processes with 

b 

^ I polarized particles are in progress. Important are setups with deuteron targets or beams 

c ■ 

m, 1^, 1^, since the deuteron serves as a unique source of information on the neutron properties 
^ , at high transferred momenta, the knowledge of which allows, e.g. to check a number of QCD 
I predictions and sum rules. Additionally, the hadron-deuteron processes can be considered 
as complementary tool in investigating short- distance phenomena and also as a source of 
information unavailable in electromagnetic reactions. Of interest is the study of nucleon 
resonances, checking non-relativistic effective models, meson-nucleon theory, NN potentials 
etc. In this line is the investigation of the deuteron break-up reaction with a fast pp pair at 
low excitation energy, proposed in and with first results reported in [^. 

One motivation for the experiment M was the possibility to investigate the off-mass shell 
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effects in A^A^ interactions. As predicted in ^, at a certain initial energy of the beam 
protons, the cross section should exhibit a deep minimum, corresponding to the node of the 
non-relativistic ^5*0 wave function of the two outgoing protons, provided the non-relativistic 
picture holds and the off-mass shell effects can be neglected. The recent data exhibits, 
however, a completely different behavior: the cross section is smoothly decreasing; there 
is no sign of a pronounced minimum. Accounting for corrections beyond the one-nucleon 
exchange mechanism improve the agreement with data, however a quantitative description 
have not achieved 0]. 

It is clear, that the non-relativistic treatment of the process becomes inadequate because 
of the high virtuality of the proton in the deuteron at the considered kinematics. More 
realistic approaches which take into account relativistic effects and the off-mass shellness 
of the interacting nucleons are desired. The Bethe-Salpeter (BS) formalism can serve as 
an appropriate approach to the problem because the off-mass shellness of the nucleons is 
an intrinsic feature of the BS equation. Moreover, the solution of the BS equation, being 
manifestly covariant, incorporate genuine relativistic effects (Lorenz boosts, negative-energy 
components etc.), hardly accessible within the Schrodinger formalism. In the present note 
we use the BS approach to analyze the data [^] on deuteron break-up with the emission of 
a fast forward pp pair 0. We pay particular attention on relativistic effects in the wave 
function of two nucleons in the continuum. The model is based on our solution of the BS 
equation for the deuteron with a realistic one-boson exchange kernel 0]. The final state 
interaction of the two protons is treated also within the BS formalism by solving the BS 
equation for the t-matrix within the one-iteration approximation P, ^j. In doing so, a big 
deal of off-mass shell effects and relativistic corrections are taken into account already within 
the spectator mechanism. 

Let us consider the process 

p +d = (pip2)(0°) +r2(180°) (1) 

at low excitation energy of the pair {E^ ~ 0— 3 MeV) and intermediate initial kinetic energies 
Tp ~ 0.6 — 2.0 GeV corresponding to the conditions at the Cooler Synchrotron COSY in 
the experiment 0, ^. In the one-nucleon exchange approximation this reaction can be 
represented by the diagram depicted in Fig. |1|, where the following notation is adopted: 
p = [Ep, p) and n = {E^, n) are the four- momenta of the incoming (beam) proton and 
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outgoing (not registered) neutron, Pf is the total four-momentum of the pp pair, which is a 
sum of the corresponding four-momenta of the detected protons, pi = {Ei, pi), p2 = {E2, P2)- 
The invariant mass of the pp-pair is M^p = P| = {2171 + ExY, where m stands for the nucleon 
mass and E^ is the excitation energy. Our calculations are performed in the laboratory 
system where the deuteron is at rest. For specific purposes, the center of mass of the pair 
will be considered as well, where all relevant quantities are superscripted with asterisks. 

A peculiarity of the processes (|l|) is that the transferred three-momentum from the initial 
proton to the second proton in the pair is rather high. Moreover, from the kinematics 
one finds that the momentum of the neutron is also high (|n| ~ 0.3 — 0.5 GeV/c), which 
implies that, since the outgoing neutron is on- mass shell, the proton inside the deuteron 
was essentially off-mass shell before the interaction. Correspondingly, it becomes clear that 
the process of A^A^ interaction in the upper part of the diagram is by far more complicate 
than the elastic interaction. For instance, let us consider a typical kinematical situation, say 
IpI = 1.22 GeV/c, e[ ~ 4°, excitation energy E^ = 3 MeV and |pi| = 0.765 GeV/c. This 
means that the neutron three-momentum is |n| ^ 0.5 GeV/c, i.e., the four-momentum of 
the off-mass shell proton was q = {Ma — En, — n). Now, if one supposes that in the upper 
vertex there was an elastic process of two on-mass shell protons, then only one kinematical 
quantity would be necessary to describe the process, e.g., at given |pi| = 0.765 GeV/c the 
scattering angle would be ~ 28° in the elastic kinematics (instead of 4° in the full reaction); 
or at given scattering angle 4*^, the momentum of the elastically scattered proton would 
correspond to |pi| = 1.21 GeV/c (|p2| = 0.334 GeV/c) instead of the detected momentum 
|pi| ~ IP2I ^ 0.765 GeV/c. This demonstrates that the A^A^ interaction in the upper part of 
the diagram has a quite complicate nature. It can be considered as consisting at least of two 
steps: (i) an inelastic process which puts the target nucleon on mass-shell, and (ii) an elastic 
interaction in the pp pair in the ^Sq final state. Since a large amount of the transferred 
energy is needed to locate the second proton on mass-shell, the relativistic corrections can 
play a crucial role here. 

The invariant cross section of the reaction (|l]) reads 

- (4?F1(M) I'^'-l'^ ' + E„) . (2) 

where X{p,d) is the flux factor, and Mfi the invariant amplitude; the statistical factor 1/2 
due to two identical particles (protons) in the final state has been already included. The 
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covariant matrix element corresponding to the diagram in Fig. |l| can be written in the form 

Mfi = u{s, n)(n - m)^d(ra) [(pa + m)^i5„(p)(pi - m)u{r, p)] , (3) 

where u{r, p) {u{s, n)) stands for the Dirac spinor of the incident proton (outgoing neutron) 
with the spin projection r (s) and 3-momentum p (n), '^£)(^isq) denote the BS amphtudes 
for the deuteron and the pp-pair in the continuum, respectively. By using the spin angular 
basis to obtain the partial decomposition of the deuteron BS amplitude and the covariant 
form for the four partial components for the ^^Q-state of the A^A^ pair |rD|, the invariant 
amplitude Mfi is 

where the contribution }C{^Sq) from the upper part of the diagram in Fig. |I] is 
ICCSo) - 



' Ej, + m 



. n|p| ^ ^ M^-2Er, ^ , , n||p| 
hi\En + m- — — - h-i \En + m + — — 

Y Ep + mJ m \ Ep + m^ 



(5) 



\l En + m 

Gs^D denote the BS vertices for the deuteron, hi, are the non- vanishing covariant partial 
components of the ^5*0 configuration in the continuum. The relation of the amplitudes 
hi {i = 1 ... 4) to the partial solutions of the BS equation in the A^A^ center of mass 
^Sq^ ,^ Sq^ ,^ ,^ Pq^ can be found, e.g., in ref. 0. In what follows we neglect the 
contribution of the extremely small component, keeping only the '^S^^ component as 

the main one and the P components as the ones providing purely relativistic corrections. 

It is easy to check that for the unpolarized particles the cross section factorizes in two 
independent parts, i.e., | — l^("'^'S'o)|^ (usi'n.)'^ + ^/^(n)^) , as it should be within 

s,r, Ai 

the spectator mechanism with ^5*0 {Lf = 0) in the final state (see also discussion in ref. [0); 
Us and ud are the BS deuteron S and D waves [|| 

In our numerical calculations we use the deuteron BS amplitude from |^. For ^5*0 state 
one should solve the inhomogeneous BS equation in the A^A^ center of mass to obtain the 
covariant amplitudes hi^^. Note that in such a way the effects of the Lorentz boost are 
automatically taken into account (see, e.g., 0]). The partial "++" components of the BS 
vertices in the A^A^ center of mass, for both the bound and the scattering states have a direct 
analogue with the corresponding non-relativistic wave functions. Moreover, at low intrinsic 
relative momenta the BS and non-relativistic wave functions basically coincide. Hence, due 
to the low excitation energy of the pp-pair, in expressing hi 3 via the partial amphtudes at 
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rest one may safely replace the "++" vertex by its non-relativistic analogue. Relativistic 
effects are then included in boosting the ^Sq^ component to the deuteron center of mass 
system and also by taking into account the contributions of P components in /ii 3. To find 
the P waves we solve the BS equation for the pair in its center of mass within the one 
iteration approximation p|, 0], with the trial function as the exact solution for the t matrix 
within a separable potential . By defining new partial components Qi {i = 1 . . .4) as 



"connected amplitudes", i.e., the partial amplitudes without the free terms, we obtain for 
the P-waves 



93{k) =ig^NN 



V3i{k,p*) - i f V3i{k,p) o p ^ • 



(6) 
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where the partial kernel V31 is V3i{k,p*) = , „ — — {\p\*Qiiy) - \k\Qo{y)} with 

|p*||k|_C/p._C/k 

as Legendre function with the argument y = (|p*P + |kp) + /i^ — A;g)/(2|p*||k|) @, In 



the first iteration the trial function gi{k) is expressed via the non-relativistic t matrix |T3 

as 

2 .c 

where the normalization of t^R corresponds to t^jiip* ,p*) = j — re* " sin 6q, where 6q 

7cm\p*\ 

denotes the experimentally known phase-shift of the elastic pp scattering in the ^5*0 state. By 
using the Sokhotsky-Weierstrass formula for the Cauchy type integrals one finally obtains 

tTTm\p*\ . * , ^ /■ r 2t/ n Jnr{p,P 



93{k) = '^'^^ 



TT 



MppV3iik,p*) 1 - '-^^^^^tr,R{p*,p*) \+mvJ dpp'Vsi{k,p) 







Ep* — Ep 



where, due to the separable stricture of the chosen t^R, further calculations of principal 
values of the relevant integrals can be carried out analytically. 

In Figs. and ^ we present results of numerical calculations of the five-fold cross section 
da / dVlidVl2d\'Pi\ and the two-fold cross section da/dQn (with fi„ as the solid angle of the 
momentum of the neutron in the center of mass of the reaction), integrated over the exci- 
tation energy in a range E^ = — 3 MeV. The calculations have been performed with our 
numerical solution for the deuteron BS amplitude (inclusion of the P components in the 
deuteron amplitude leads to negligibly small corrections, therefore here we do not discus 
these contributions). The dotted curves are results of non-relativistic calculations, while the 
dashed curves include pure Lorentz boos effects, i.e., relativistic calculations with including 
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the ^Sq~^ component only. It is clearly seen that the boost corrections are fairly visible: 
They cause a shift of the minimum of the cross section. The agreement at low initial ener- 
gies becomes better, however, the cross section is still too small at large values of Tp. Fig. |^ 
reveals that an account for only " + +" components is not sufficient to describe data. An 
excellent description is achieved by taking into account all the relativistic effects, including 
the contribution of the negative-energy P waves. Note, that as demonstrated in refs. |^ 
in reactions of pd and in threshold-near ed disintegration, the inclusion of P waves exactly 
recovers the non-relativistic calculation with taking into account the NN pair production 
effects. Hence, in our case this is a hint that covariant calculations within the relativistic 
spectator mechanism contains already some contributions beyond the one-nucleon exchange 
mechanism. 

In summary we analyze the recent data ||^ of the reactions pd{p, n) —>■ {pp)n within 
a covariant approach based on the Bethe-Salpeter formalism. Relativistic effects (Lorentz 
boost, negative-energy P components) are important and responsible for the smooth decline 
of the cross section, in contrast to predictions of non-relativistic models. 
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FIG. 1: Kinematics of the process (|l|). 
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FIG. 2: Five-fold cross section as a function of the kinetic energy Tp of the incident proton. The 
dotted curve corresponds to a non-relativistic calculation, i.e., to the case where only the "+- 1-" 
components in the state are taken into account and any Lorentz boost effects ignored. The 
dashed curve depicts results of calculations with all relativistic effects in "+- 1-" components taken 
into account. The solid curve is for the results of a complete calculation with taking into account 
all the relativistic effects including the contribution of P waves in the wave function of the pp pair. 
The two protons are supposed to be detected in the strictly forward direction, i.e., 9\ = O2 = ■ 
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FIG. 3: Differential cross section in the center of mass integrated over the excitation energy Ex as 
a function of the kinetic energy Tp of the incident proton. Notation is as in Fig. |2[ Experimental 
data are from 



